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1.1  Objective 
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Before  deploying  e  new  veep on  system,  the  US  Army  must  determine  how  to 
support  the  system  In  the  most  cost  effective  manner.  Currently,  sophisticated 
multi-echelon  models  compute  stockage  quantities  for  spares  and  repair  parts 
which  minimize  total  inventory  investment  while  achieving  a  target  level 
of  operational  availability.  The  maintenance  policies  to  be  followed  are  input 
to  the  stockage  models. 

The  "Optimum  Allocation  of  Test  Equipment /Manpower  Evaluated  Against 
Logistics"  (OATMEAL)  model  determines  optimum  maintenance  as  well  as  stochage 
policies.  Specifically,  it  determines  at  which  echelon  each  maintenance  function 
should  he  performed,  or  whether  the  maintenance  function  should  be  eliminated;  i 
it  does  repair  vs  throwaway  analysis  as  well  as  level  of  repair  analysis. 

The  costa  minimized  include  those  which  depend  on  the  range,  number  and 
placement  of  test  equipments,  those  which  depend  on  whether  repair  is  done  anu 
where  it  is  performed,  and  those  which  depend  on  the  range  and  dollar  value  of 
repair  parts  and  assemblies  stocked  at  eech  echelon.  The  model  calculates  inven¬ 
tory  levels  and  maintenance  policies  which  will  achieve  the  operational  avail¬ 
ability  target  Input  by  the  user.  It  will  not  choose  maintenance  policies  that 
preclude  achievement  of  the  target. 

The  objective  is  to  achieve  the  target  operational  availability  at  minimum 
life  cycle  cost. 

1.2  Conceptualization  of  Maintenance 

The  model  looks  at  three  levels  of  indenture  within  a  weapon  system;  com¬ 
ponents,  modules,  and  piece  parts.  The  indenture  breakdown  of  a  home  audio  system 
is  depicted  in  Figure  1.  A  failure  mode  is  defined  as  a  system  failure  due  to  « 
specific  module  in  a  specific  component.  All  failure  rates  are  input  by  failuie 
mode,  and  maintenance  allocation  decisions  are  also  made  by  failure  mode.  OATMEAI. 
considers  four  echelons  of  maintenance:  (1)  organizational,  (2)  direct  support  (1>?) , 
(3)  general  support  (GS) ,  and  (4)  depot . 

Consider  maintenance  on  the  audio  system  shown  in  Figure  1.  In  every  case  of 
turntable  failure,  the  audio  system  is  repaired  at  ORG.  If  the  arm  mechanism  fails, 
the  ORG  replaces  the  turntable  and  returns  the  unserviceable  turntable  to  the  USi1. 


£ 

I 


1 1*.  i  •*,  »*.  ’* 

^  ■  ■»  t ..  +  i  ll 


•  •  »  '  »  ^  ’  v*  V  • 

h  '.s,  t  m  .*  a 


At  the  DSU  a  new  arm  mechanism  la  put  in,  end  the  old  arm  mechanism  Is  thrown 
away.  Tf  the  motor  falls*  the  unserviceable  turntable  will  be  evacuated  all  the 
way  to  Gfitl  for  repair.  After  the  motor  la  replaced  at  GSl>,  the  unserviceable  motor 
La  shipped  to  depot  for  repair.  If  the  needle  fella*  all  work  is  done  at  ORG,  and 
the  needle  is  thrown  away. 

The  user  any  have  originally  daalgnatad  the  needle  aa  a  consumable,  but 
the  arm  mechanism  as  pctantially  rspalrablc*  whlla  tha  model  concluded  it  was 
not  cost  effective  to  repair  the  arm. 

It  sometimes  happens  that  the  aaae  module  appears  in  two  (or  more)  different 
components.  Failure  rates  would  be  luput  separately  for  each  module  applica¬ 
tion,  but  the  commonality  would  ba  racognlsed  in  computing  atockage  and  logisti¬ 
cal  costs. 

Since  detailed  pert  data  is  not  generally  available  in  early  development, 
the  pieceparts  are  considered  in  an  aggregate  manner.  It  is  also  possible  to 
represent  a  group  of  modules,  or  even  components  by  an  average  module  or  com¬ 
ponent,  specifying  how  many  distinct  modules  or  components  this  average 
represents. 

1 • 3  Test  Equipment  and  Repair  Skills 

The  need  to  use  test  equipment  which  can  be  quite  expensive  makes  the 
level  of  reDair  analysis  mathematically  challenging.  A  given  repair  action 
may  require  more  than  ->ne  piece  of  equipment,  and  many  different  actions  may 
have  a  requirement  for  the  same  piece  of  equipment .  It  is  assumed  that  any 
piece  of  equipment  required  for  fault  diagnosis  will  also  be  required  for 
repair,  so  that  diagi.osis/repalr  can  be  considered  one  procedure.  This  is 
reasonable  because  usually  a  repair  is  not  considered  complete  unless  the 
equipment  used  for  diagnosis  has  performed  a  functional  check  to  verify  the 
success  of  the  repair. 

Test  equipment  is  labelled  as  either  common  or  peculiar  at  each  echelon. 

If  peculiar,  only  Integer  quantities  cat.  be  placed  at  a  repair  facility.  The 
whole  cost  of  the  equipment  must  be  considered,  even  though  the  equipment  is 
used  at  tha  facility  only  a  fraction  of  the  time.  If  the  test  equipment  is 
common,  it  Is  assumed  that  only  use  must  be  paid  fcr;  if  the  end  item  needs 
l/4  of  the  throughput  of  a  common  piece  of  test  equipment  at  the  facility,  it 
bears  1/4  of  the  total  cost  for  the  equipment. 

Test  equipment  may  be  needed  for  three  different  types  of  repair  actions: 

a.  Repair  weapon  system  when  it  fails  due  to  the  failure  of  a  specific 


component. 
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b.  Repair  a  component  whan  it  fails  due  to  the  failure  of  a  specific  modulo. 

c.  Repair  a  nodule  vhenever  it  faila. 

Repair  Skills.  Coat  of  repair  including  labor  costs  Is  Input  tc  the  model . 
However,  sonatinas  repair  require*  apodal  skills  which  would  be  new  to  an  echelon. 
The  user  nay  identify  special  repair  skills,  so  that  the  training  costs  and  higher 
salaries  are  properly  calculated.  If  a  special  skill  is  needed  at  DSU,  hut  the 
skill  will  be  required  for  only  1/4  of  e  nan  year  per  year,  full  training  costs 
for  one  person  per  DSU  must  still  bs  incurred.  The  model  accounts  for  tills  just 
as  it  properly  accounts  for  paculiar  test  equipment.  In  both  rases  there  is  an 
incentive  to  allocate  maintenance  ao  that  skills  and  equipments  are  not  needed 
at  echalons  where  they  will  bo  grossly  underutilised. 

1.4  Inputs  end  Outputs 

Input  consists  of  ths  type  of  data  neceasary  to  run  a  supply  model  such 
os  SESAME  [10],  plus  additional  cost  and  maintenance  related  data.  The  user 
can  also  provide  Inputs  which  exclude  certsin  types  of  solutions  as  not  feas¬ 
ible  or  raalistic;  for  example,  he  can  specify  that  a  certsin  test  equipment 
cannot  be  placed  below  the  GSU.  The  maintenance  data  consists  primarily  of 
test  equipment  end  special  skill  requirements  for  each  maintenance  action. 

There  can  be  up  to  four  eupply  and  maintenance  echelons  Allowed.  It  is 
assumed  supply  and  maintenance  functions  are  colocated,  but  an  intermediate 
echelon  can  be  specified  to  have  only  a  maintenance  function. 

Replacement  task  distributions  end  maintenance  task  distributions  [11] 
are  outputs  of  OATMEAL  whereas  they  would  ba  input  to  supply  models.  An  item's 
replacement  task  distribution  gives  the  percent  of  removals  which  occur  at  each 
echelon,  while  the  maintenance  task  distribution  specifies  the  percent  of  removed 
items  which  can  be  repaired  at  each  echelon. 

Other  outputs  of  OATMEAL  are  the  numbers  and  locations  of  test  equipment 
and  special  personnel,  end  the  stockage  quantities  associated  with  the  main¬ 
tenance  policy  chosen. 

1.5  Solution  Approach 

The  heart  of  OATMEAL  is  a  mixed  integer  program  (MIP)  mode  .  Consider  the 
approach  to  solution  to  consist  of  four  stages  -  preprocessing,  formulation, 
optimization,  and  evaluatlor .  The  solution  flow  is  depicted  in  Figure  2. 

Preprocessor .  This  program  [8]  accepts  input  from  the  user  in  a  format 


designed  to  be  as  convenient  as  possible.  Default  values  for  inputs  are  pro¬ 
vided  whenever  it  makes  sense  to  do  so,  and  thera  are  edit  checks.  The  output 
is  a  parsimonious  set  of  parameter  and  variable  values  necessary  to  describe 
the  problem  mathematically. 

Formulator.  This  program  formulates  the  problem  as  a  MIP  model  in  the 
specific  format  for  the  commercial  software  package  being  used  (see  below). 
Ucwever,  to  do  this.  It  must  first  generate  stockage  (and  cost  Implication 
thereof)  for  all  candidate  maintenance  policies.  The  Formulator  uses  the  SESAME 
mult ’.echelon  stockage  modal  [5]  as  a  subroutine  to  calculate  optimally  quantities 
of  components  and  modules  for  a  particular  maintenance  policy  vector.  By  policy 
vector  wo  mean  an  assignment,  for  each  failure  mode,  of  where  to  replace  and  where 
to  repair  the  respective  components  and  modules.  The  associated  stockage  costs 
will  then  be  accessible  to  the  MIP  cost  function  formulation. 

Optimizer;  The  "MIP'*  refers  to  the  commercial  software  package  which 
accomplishes  the  optimisation.  Packages  for  accomplishing  MIP  are  available,  each 
of  which  required  many  man  years  to  develop.  The  specific  package  now  being  used 
is  APEX  III,  (1),  developed  by  Control  Pate  Corporation,  but  others,  such  as  IBM's 
KPSX,  can  be  substituted  with  little  change. 

MIP  "optimizes"  ••  finds  the  maintenance  policy  vector  which 
minimizes  costs,  including  stockage  costs  and  backorder  penalty  costs  *  by  using 
an  efficient,  heuristic  search  so  that  all  the  myriad  combinations  do  not  have 
tc  be  evaluated. 

Eve luator.  As  the  name,  implies,  this  program  evaluates;  it  accepts  as  input 
all  the  data  about  the  problem  being  run  and  also  a  maintenance  policy  vector 
(the  replace-repair  actions  by  failure  mode)  to  be  assessed  In  terms  of  cost  and 
performance.  This  policy  vector  may  have  come  from  the  optimizer  or  been  proposed 
bv  the  user  (rote  alternative  policy  block  in  Figure  2) .  The  evaluator  uses 
SESAME  subroutine  to  determine  for  each  component  end  module  the  optimum  stockage 
quantities.  The  evaluator  will  determine  the  operational  availability  performance 
based  on  the  computed  stockage  quantities  and  the  selected  maintenance  policy 
vector;  it  also  will  compute  all  relevant  costs  (see  Chapter  2  for  cost  elements 
and  expressions).  Finally  the  evaluator  will  convert  the  failure  mode  policy 
vector  to  component  and  module  replacement  task  distributions  and  maintenance 
task  distributions  (Section  A. 8.4). 

OATMEAL,  by  incorporating  SESAME  model  subroutines,  simultaneously  optimizes 
maintenance  and  supply..  If  SESAME  or  a  comparable  stockage  optimizer  is 
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not  included  in  a  repair  level  analysis  tiodel,  eetlaatee  of  the  impact  of  different 
maintenance  policies  on  atockage  coats  cannot  bt  correctly  estimated.  Therefore, 
this  nay  result  ir  the  choice  of  maintenance  policies  which  are  not  as  cost  effec¬ 
tive  as  expected.  This  phenomenon  will  be  most  pronounced  when  the  maintenance 
policies  have  an  adveraa  impact  on  operational  availability,  requiring  a  very 
large  investment  in  inventory  to  compensate. 

1.6  Docunantation  of  OATMEAL. 

The  User's  Guide  (8)  explains  la,  detail  the  input  required  from  the  users, 
default  parameters,  the  interpretation  of  output,  and  different  approaches  t 
using  the  OATMKAL  model  for  maximum  benefit.  It  also  documents  the  transforma¬ 
tions  made  by  the  pre-processor  to  develop  the  inputs  required  by  the  MIP  foi mu¬ 
lct  or. 

This  report  is  Intended  for  a  range  of  readers  with  different  objectives. 

This  Chapter,  Chapter  £ I  and  Sections  3. 1-3. 3  of  Chapter  III  do  not  require  a 
sophisticated  mathematical  background,  and  are  Intended  to  give  analytically 
Inclined  readers  e  good  understanding  of  the  capabilities  of  the  model.  Chapter 
II  describee  how  the  model  evaluates  costa  end  operational  svailebility.  Chapter 
HI  discusses  the  elements  of  the  optimisation  process,  i.e.  the  MIP  formulation, 
selection  of  e  backorder  penalty  parameter,  end  evaluation  of  the  cost-availability 
performance  for  maintenance  policy  vectors. 
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MAINTENANCE  CONCEPTS 


FIGURE  1 


CHAPTER  IX 


EVALUATION  OF  COST  AMD  AVAILABILITY 

2.1  Introduction 

This  chapter  will  give  the  reader  the  scope  of  the  costs  considered,  the 
actual  cost  equations,  and  an  idoa  of  how  performance  is  evaluated  through  end 
item  availability.  A  more  mathematical  treatment  of  the  functional  details  of 
OATMEAL  la  left  to  Chapter  III. 

2.2  Use  of  Present  Value 

Because  OATMEAL' s  objective  la  to  minimize  life  cycle  costa,  the  costs  con¬ 
sidered  are  a  mixture  of  one-time  costa  and  recurring  annual  costs.  To  make 
these  two  types  of  costs  commensurable,  the  present  value  approach  as  recom¬ 
mended  in  DoDI  7041.2  [12]  la  used. 

As  an  example,  given  that: 

a.  100  requialtiona  are  processed  per  year. 

b.  Coat  per  requisition  processed  is  $10. 

c.  Expected  lifetime  for  the  system  Is  11  years. 

d.  Discount  rate  for  present  value  analysis  is  102. 

Then,  assuming  costs  are  Incurred  at  midyear, 

11  ,  t-1/2 

Requisition  Costs  (Pressnt  Value)*  E  C$10)  x  (100)  x  (■  y  77:) 

t-1  1,10 

This  may  be  rewritten  as: 

11  .  t-1/2 

Requisition  Costa  (Present  Value)  -  ($10)  x  (100)  x  -  (T~yn) 

t-1  1,10 

-  Annual  Requisition  Costs  x  PVFAC 

where  PVFAC  -  present  value  adjuetment  factor.  Note  that  PVFAC  depends  only  on 
the  expected  life  of  the  system  and  the  discount  rate,  and  can  be  used  to  con¬ 
vert  other  annual  coats  to  the  present  value  of  expected  lifetime  costs.  While 
PVFAC  may  be  calculated  as  shown,  it  is  input  to  OATMEAL  from  the  pre-processor, 
which  takes  the  correct  value  from  official  DoD  tables. 

2 . 3  Full  Deployment  Assumptions 

Full  deployment  f.s  assumed  in  year  1.  This  means  cost  estimates  produced 
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are  not  true  life  cycle  costs,  but  are  useful  for  ranking  policy  alternatives. 

This  approach  was  taken  to  simplify  data  requirements  and  processing  with  the 
expectation  that  this  would  not  unduly  bias  choice  of  one  alternative  over 
another.  It  does  tend  to  exaggerate  the  impact  of  costs  which  will  phase  in 
over  tine,  so  some  refinements  may  ultimately  be  necessary,  e.g.,  a  factor 
applied  to  PVFAC,  Note  that  not  only  do  annual  costa  such  as  requisition  pro¬ 
cessing  build  up  to  full  deployment  levels,  but  even  "onetime''  costs  such  as 
purchase  price  for  test  equipment,  are  not  really  all  incurred  at  one  time; 
test  equipment,  for  example,  need  only  be  deployed  as  the  weapon  system  is 
introduced  over  time  to  additional  fighting  units. 

2 • 4  Logistical  Support  Cost  Components 

Costs  are  computed  for  each  individual  component  and  module  using  the  equa¬ 
tions  described.  The  variables  underlined  in  the  equations  are  passed  directly 
from  the  preprocessor.  The  other  variables  are  computed  or  modified  by  OATMEAL 
and  depend  on  the  maintenance  policies  chosen. 

2.4.1  Throw  Out  Costs 

Throw  out  costs  represent  the  annual  value  of  the  components  and  modules 
which  wash  out  and  must  therefore  be  replaced  by  nev  procurement.  As  with  all  the 
ether  annual  costs,  PVFAC  is  used  to  convert  the  estimate  of  annual  costs  to  the 
present  value  of  expected  lifetime  throw  away  costs.  It  is  assumed  that  the 
administrative  costa  of  making  procurements  will  not  significantly  vary  among 
alternative  maintenance  policies  and  need  not  be  considered. 

Mathematically,  the  equation  for  throw  out  costs  for  any  given  component 
or  module  is : 

THROW  OUT  COST  -  FVFAC  x  (Annual  Removals)  x  (%  Washout)  x  (Unit  Price) 

Annual  removals  reflect  failure  rates  plus  "false"  removals.  A  "false" 
removal  occurs  when  a  component  or  module  which  is  perfectly  good  is  removed  due 
to  an  error  or  ambiguity  in  diagnosis.  Currently,  the  model  user  inputs  a  failure 
rate  for  each  failure  mode,  but  a  single  false  removal  rate  which  applies  to  all 
items.  If  the  false  removal  rate  is  10  percent  and  a  component  or  module  is  ex¬ 
pected  to  fall  100  times  a  year,  total  removals  are  estimated  as  (100  +  100  x  10%) 
or  110. 

Maintenance  policies  impact  removal  rates  in  that  if  a  component  is  thrown 
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out  rather  than  repaired,  this  will  eliminate  demand  for  modules  used  to  repair 
that  component.  Throw  out  coats  are  computed  for  the  component,  since  each  com¬ 
ponent  thrown  out  must  be  replaced  by  procurement,  but  no  throw  out  costs  are 
computed  for  its  modules. 

At.  Inherent  washout  rate  is  input  for  each  component  and  module,  reflecting 
the  percent  of  removals  which  cannot  be  fixed  regardless  of  the  maintenance 
policies  chosen.  If  a  throw  out  policy  is  selected  by  the  model,  or  input  for 
evaluation,  the  washout  rate  becomes  100  percent. 


2.4.2  Repair  Costs:  Conson  Labor  and  Manuals 


Repair  costs  include  labor,  parta,  and  the  need  to  develop  repair 
manuals.  Two  types  of  labor  are  considered,  common  and  special  skills.  The 
reason  for  treating  special  skills  separately  is  that  just  as  peculiar  test 
equipment,  they  may  not  be  fully  utilized;  e.g.,  if  repair  requiring  a  special 
skill  is  done  at  ORG,  the  full  cost  of  putting  the  specially  trained  person  at 
ORG  is  incurred  even  though  he  may  need  to  use  his  special  skill  only  a  small 
portion  of  the  time.  Special  skills  costs  are  discussed  further  in  Section  2.5 
while  parts  costs  are  discussed  in  Section  2.4.8. 

While  special  skill  requirements  are  treated  in  detail,  common  labor  costs 
are  treated  somewhat  approximately  in  order  to  simplify  the  data  requirements 
of  the  model.  The  user  Inputs  the  average  hours  of  labor  to  repair  each  component 
and  each  module.  He  does  not  enter  the  hours  to  repair  the  end  item  itself,  or 
relate  the  number  of  hours  to  repair  a  component  to  the  failure  mode.  Thus, 
common  labor  repair  costs  are  computed  per  component  or  per  module  repair  action. 

Because  of  differences  in  pay  scales  and  working  hours  by  echelon,  the  same 

job  will  incur  a  different  cost  depending  on  the  echelon  at  which  it  is  done. 

*» 

The  pre-processor  takes  all  this  information  into  account  and  inputs  to  OATMEAL 
the  cost  per  repair  action  by  echelon.  The  number  of  repair  actions  per  echelon 
is  then  computed  by  the  EVALUATOR  from  the  annual  removals  and  the  maintenance 
task  distribution. 


MT^k  as  the  porcent  of  all  removals  (no  matter  where  removed) 

repaired  at^Khalon  k,  where  k  takes  on  the  values  1,  2,  3  and  4  to  refer 

respectively  to  the  ORG,  DSU,  GSU  and  depot  repair  echelons  (e.g.  MTD^  is  percent 

of  all  removals  repaired  at  GSU),  The  oum  of  MTD.  +  MTD„  +  MTD.,  +  MTD,  plus  the 

a  l  3  4 

washout  rate  will  always  equal  100  percent. 


Repair  Costfc  la  defined  as  ths  cost  par  repelr  action  if  dons  at  echelon  k. 
Mathematic ally ,  for  each  component  and  module 

COMMON  LABOR  COSTS  -  (PVFAC)  x  (Annual  Removals)  x  £ (MTD^)  x  (Repair  Costk) 

Preparation  of  manuals  is  viewed  as  a  composite  of  documenting  each  component 
and  each  module  repair  action.  The  pre-processor  estimates  these  costs  based  on 
the  number  of  pages  required  to  document  each  type  of  repair  action  as  input  by 
the  user.  The  EVALUATOR  adds  manual  preparation  cost  to  component  and  module 
repair  cost  unlass  the  policy  is  to  do  no  repair  (throwaway)  of  the  component 
or  module. 

2.4.3  Transportation  Costa 

Transportation  cost  evaluations  are  straightforward  once  we  clarify 
the  distinctions  between  demands,  removals,  and  retrograde.  He  do  this  with  an 
example  chosen  for  its  pedogoglcal  value  rather  than  its  realism.  Suppose  we 
have  a  simple  supply  system  consisting  of  organisations,  direct  support  units, 
and  a  depot.  For  a  particular  module  (denoted  PMOD)  there  are: 

30  ORG  echelon  PMOD  removals  per  year 
10  ORG  echelon  PMOD  repairs  per  year 
3  ORG  echelon  PMOD  washouts  per  year 
50  DSU  echelon  PMOD  removals  per  year 
40  DSU  echelon  PMOD  repairs  per  year 
5  DSU  echelon  PMOD  washouts  per  year 

The  ORGs  will  demand  20  PMODs  a  year  from  the  DSUs  (30  removals  minus  10  repairs) . 

The  DSU's  need  to  supply  these  20  and  replace  the  50  PMODs  removed  qg  DSU.  Since 

the  DSUs  are  able  to  repair  40  a  year,  their  net  shortfall  is  (50  +  20  -  40)  or 

30.  Thus  at  DSU  level  there  are  50  removals  a  year  and  30  demands  on  the  depot. 

The  ORGs  are  retrograding  (30  -  10  -  3)  or  17  PMODs  to  the  DSUs  for  repair. 

The  DSUs  see  these  17  unserviceable  PMOD's  plus  50  removals  at  DSU,  and  then 

repair  40  and  washout  5.  So  DSU  retrograde  to  the  depot  is  (17  4*  50  “  5) 

or  22.  Currently,  two  assumptions  are  made  which  relate  to  transportat^bn  costs. 

It  is  assumed  only  repairable  items  are  retrograded.  Thus,  in  our  example  we 

1 

assumed  the  ORGs  were  knowledgeable  enough  to  Identify  three  PMODs  as  washout 
and  not  ship  them  to  the  DSU  with  the  17  PMODs  retrograded.  It  is  also  assumed 
retrograde  cost  rates  are  equal  to  forward  transportation  rates.  For  more 
details  or  retrograde  calculations  see  Appendix  A  Section  A. 8. 3. 
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The  reader  can  sea,  from  the  example  given,  that  the  maintenance  policies, 
by  affecting  where  repair  la  done,  l  both  demand  and  retrograde.  Let 

Anrual  Demand^  denote  the  total  demands  placed  by  echelon  k  sites  on  their 
suppliers  and  Transport  Rate^  denote  the  cost  per  pound  of  moving  material 

and  Retro  Rate^ 

refer  to  movement  from  echelon  k  to  the  maintenance  units  providing  support. 

Then, 

TRANSPORTATION  COST  - 

PVFAC  x  (Item  Weight)  x  I  (Annual  Demand^)  x  Transport  Rate^) 

+  PVFAC  x  (Item  Weight)  x  E  (Annual  Retrograde^)  x  (Retro  Rate^) 

k 

2.4.4  Requisition  Coats 

Annual  Requisitions^  is  defined  as  the  number  of  requisitions  for 
a  given  component  or  module  placed  annually  by  echelon  k  sites  on  their  suppliers. 
One- for- one  ordering  is  assumed  for  the  moat  part,  i.e.  each  item  demanded  results 
in  an  additional  requisition.  However,  this  is  qualified  to  the  extent  that  it 
is  assumed  no  single  site  will  requisition  the  same  item  more  than  12  times  a 
year.  Thus,  for  low  demand  items,  number  of  requisitions  equals  number  of  demands, 
but  not  for  high  demand  items,  for  which  Annual  Requisitions^  may  not  exceed  12. 
Mathematically, 

REQUISITION  COSTS  ■  PVFAC  x  l  (Annual  Requisitions.  )  x  (Requisition  Cost  Parameter) 

k 

2.4.5  Stockage  Costs 

Stockage  coats  include  the  one-time  costs  of  buying  stock  to  fill 
the  pipeline,  and  an  annual  holding  cost  to  cover  stockage  and  losses  or  pilferage 
of  inventory.  One  issue  is:  what  is  the  impact  on  cost  of  engineering  redesign 
whereby  better  performing  components  or  modules  replace  those  previously  used? 

Will  the  new  item  simply  be  bought  when  it  is  time  to  replenish  washouts  of  the 
old  item,  or  will  all  stocks  of  the  old  item  be  excessed  at  time  of  engineering 
redesign,  increasing  cost?  A  second  issue  is:  Can  pipelines  be  drawn  down 
gradually  so  that  if,  for  example,  phase  out  begins  in  1994,  procurements  to 
replace  washouts  after  1994  are  avoided  by  using  stocks  in  the  pipeline,  saving 
money? 


from  supplier  to  echelon  k.  Analogously,  Annual  Retrograde^ 
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To  simplify  data  requirements  the  laauaa  raised  are  avoided  by  not  Including 
either  excess  costs ,  or  swinge  from  pipeline  reduc  tions,  letting  HCFAC  be  the 
annual  cost  of  storage  and  pilferage,  as  a  percent:  of  unit  price: 

STOCKAGE  COST  -  (Quantity  Stocked)  x  (Unit  Price)  + 

PVFAC  x  (Quantity  Stocked)  x  (Unit  Price)  x  HCFAC 

2.4.6  Bln  Costs 

Bln  costa  are  those  management  and  holding  costa  which  vary  as  a 
function  of  the  range  uf  items  stocked  rather  than  the  dollar  value  or  quantities 
stocked.  The  bln  cost  parameter  is  the  cost  per  NSN  per  atockage  location  per 
year  so  that  for  each  component  or  module: 

BIN  COST  •  PVFAC  x  (Number  of  Stocking  Locations)  x  (Bln  Cost  Parameter) 

Bln  costs  and  stockage  costs  are  the  two  cost  components  that  depend  on  the 
answers  found  by  the  SESAME  subroutines  of  OATMEAL  as  to  how  much  and  where  to 
stock  components  and  modules. 

2.4.7  Catalog  Costs 

Each  module  and  component  is  coded  as  to  whether  It  is  a  new  item  or 
not.  If  It  is  new,  a  one  time  catalog  Introduction  cost  is  incurred;  additionally 
a  recurring  maintenance  cost  is  assessed. 

Therefore, 

CATALOG  COSTS  ■  (Total  New  Items)  x  (Item  Introduction  Cost)  + 

PVFAC  x  (Total  New  Items)  x  (Item  Maintenance  Cost) 

2.4.8  Parts  Costs 

For  each  module  an  average  part  is  created  to  represent  all  parts 

used  in  fixing  the  module.  The  unit  price  to  be  used  for  the  average  part  is 

calculated  by  the  pre-processor  based  on  the  average  value  of  parts  used  per 

repair  action.  Demand  is  estimated  as  the  total  demand  for  the  module,  less  module 

* 

washouts  divided  by  the  number  of  parts  the  average  represents.  The  echelons 
at  which  demand  arises  are  inferred  from  the  echelons  at  which  tne  module  is 
repaired. 

* 

The  number  of  parts  is  actually  baaed  on  new  parts  only.  This  means  demand  for 
old  parts  is  attributed  to  new  parts  so  that  added  stockage  costs  for  old  parts 
are  reflected  an  added  cost  for  new  parts. 
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Once  the  average  part  la  described,  SESAME  subroutines  are  used  to  calculate 
atockage  and  all  logistics  costs  are  calculated  Including  throwaway  (washout  is 
always  lOo  percent),  requisitioning  cost,  stockage  cc  it  and  bln  cost.  These  coats 
are  multiplied  by  the  number  of  parts  the  average  represents.  The  number  of  new 
parts  Introduced  contributes  to  cataloging  costs. 

2.4.9  Backorder  Ccfti 

Component  backorders  degrade  operational  availability,  while  module 
backorders  Increase  the  number  of  components  in  the  repair  pipeline.  Ideally, 
these  Increased  pipelines  would  be  accounted  for  in  computing  component  stockage. 
Neither  SESAME  nor  the  OATMEAL  model  currently  does  this .  Therefore ,  the  back¬ 
orders  themselves  are  costed  out.  If  module  backorders  will  Increase  the  number 
of  components  in  the  pipeline  by  n,  and  each  component  costs  UP,  module  backorders 

A 

are  costed  out  as  (n)  x  CUP)  < 

By  analogous  reasoning,  parts  backorders  are  costed  out  in  terms  of  the 
Increase  in  the  expected  number  of  modules  in  the  repair  pipeline. 

2.5  Teat  Equipment  and  Special  Manpower  Coats 

For  each  type  of  test  equipment  and  each  different  kind  of  specially  tra  .ed 
repairman,  the  EVALUATOR  calculates  the  total  requirement  and  multiplies  by  the 
cost  per  equipment  or  per  repairman  which  is  input.  The  pre-processor  bases 
test  equipment  cost  on  purchase  price,  installation  cost,  and  an  annual  main¬ 
tenance  cost  expressed  as  a  percent  of  purchase  price.  It  also  checks  to  see  if 
the  test  equipment  life  is  less  than  the  weapon  system's  life,  in  which  case  it 
includes  a  cost  to  reflect  the  need  to  eventually  replace  the  test  equipment. 

The  pre-processor  bases  repairman  costs  on  salary  and  training  cost;  annual 
training  cost  is  the  cost  of  training  divided  by  the  average  length  of  time  a 
repairman  stays  in  his  position.  Repairman  cost  may  vary  by  echelon 

OATMEAL  calculates  workload  on  equipment  or  special  repairmen  in  detail. 
Workload  factors  are  input  to  the  EVALUATOR  by  failure  mode  for  end  itati  and 
component  repair;  l.e.  workload  to  repair  the  end  item  may  depend  on  which 
component  failed,  and  workload  to  fix  the  component  may  depend  on  which  module 
failed.  Module  repair  workloads  do  not  depend  on  which  component  the  module  came 
from,  nor  does  the  user,  in  developing  input  to  the  pre-processor,  necessarily 


Actually,  to  be  consistent  with  calculation  of  stockage  cost,  we  use 
(n)  x  (UP)  +  (PVFAC)  x  (n)  x  (UP)  x  (HCFAC) . 
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have  to  use  the  capability  to  aak*  component  and  end  Iten  repair  workloads  depend 
on  failure  node. 

Equlpaent  and  repairmen  requlr  'nts  are  calculated  by  site,  reflecting 
the  workload  factors,  pqulpmants  r„  ^rted  by  the  site,  and  the  maintenance 
policies.  If  the  equipment  or  repairman  will  be  at  the  site  just  to  support 
the  single  weapon  system,  l.e.  It  lu  peculiar  to  the  weapon  system,  requirements 
are  rounded  up:  2.2  becaams  3  and  so  on. 

2 . 6  Evaluation  or  Operstlonal  Availability 

Operational  availability  la  estimated  as: 

n.  MCTBF _ 

"  HGT1F  +  MIR  4*  NIT  +  MLDT 

where 

OA  -  operational  availability  of  the  weapon  system 

MCT1T  -  mean  calendar  time  between  failures 

MTR  -  mean  time  to  repair  the  weapon  system  If  all  resources  are 
available 

MTT  -  mean  transportation  time 

MLDT  -  mean  logistics  down  time 

MCTBF  and  MTR  are  Inputs  and  relate  to  the  perfoneance  of  the  weapon  system 
Independent  of  what  support  It  receives. 

MTT  Is  a  function  of  the  repair  level  decisions.  If  the  ayetem  Is  always  repaired 
at  user  level  with  user  personnel  and  equipment,  MTT  equals  0.  Otherwise,  MTT 
covers  the  time  for  the  upper  echelon  personnel  to  get  to  the  user,  or  for  the 
system  to  be  moved  to  the  repair  site  and  back. 

MLDT  la  the  mean  time  to  get  an  essential  component  from  the  supply  system 
when  needed  for  weapon  system  repair.  It  depends  on  the  repair  level  analysis 
in  that  for  any  given  set  of  maintenance  policies,  there  is  an  associated  set  of 
stockage  levels,  and  the  MLDT  is  a  function  of  these  supply  levels  and  the 
maintenance  policies.  The  SESAME  stockage  model  calculates  MLDT. 
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OPTIMIZATION 

3*1  Introduction 

This  chapter  discuaaea  the  mathematics  and  procaduraa  for  problem  formula- 
tlon  and  optimisation. 

One  difficulty  encountered  relate*  to  the  need  to  use  SESAME  in  problem 
formulation.  To  compute  stockage,  SESAME  muat  first  find  the  value  of  a 
Lagranglan  (called  CURPAR  in  the  SESAME  literature)  by  which  stockage  Is  related 
to  the  operational  availability  target.  Unfortunately,  SESAME  cannot  find  the 
Lagranglan  without  first  knowing  what  maintenance  allocation  decisions  will  be 
made. 

Thus,  there  la  circularity:  optimimum  maintenance  allocation  decisions 
depend  on  atockaga  coats,  while  the  computation  of  stockage  quantities  require 
that  CURPAR  be  known,  and  CURPAR  depends  on  what  the  maintenance  allocation  de¬ 
cisions  are.  In  Sections  3.5  and  3.6  wa  describe  the  circularity  and  how  we 
circumvent  it  in  detail.  Until  then  the  reader  is  asked  to  accept  that  SESAME 
can  compute  truckage  in  the  problem  formulation  stage. 

3.2  Problem  Formulation  Procedures 

Aa  Jtated  in  Chanter  I,  OATMEAL  considers  three  levels  of  indenture  and 
four  echelons  of  maintenar ca  (pluu  throw-away) .  This  structure  would  generate 
5^  possible  maintenance  policies  for  each  failure  mode.  However ,  OATMEAL  re¬ 
quires  that  a  modi' la  be  repaired  at  the  same  or  higher  echelon  as  the  component, 
and  the  component  repaired  at  the  ernes  or  higher  echelon  as  the  end  item. 

Further,  aud  item  repair  ia  limited  to  ORG  or  DS,  resulting  in  25  possible 
maintenance  teak  allocations  for  * ich  failure  mode.  These  policies  can  be  ex¬ 
pressed  as  a  triplet  in  which  the  first  variable  states  the  echelon  for  end  item 
repair,  the  second,  component  repair,  and  the  third,  module  repair  (see  Table  3.1). 
In  particular  cascr,  soma  of  the  25  may  be  excluded  by  the  user.  Although  the 
optimisation  examines  25  maintenance  policies  for  each  failure  mode,  the  formu- 
lator  does  not.  Instead,  the  costa  for  these  policies  are  built  up  from  the 
costs  found  for  components  and  modules  conaldsred  individually.  Thus,  in  computing 
stockage,  the  formulrtor  applies  SESAME  to  each  of  nine  alternatives  for  each 
component  end  to  each  of  15  alternatives  for  each  module. 
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COMPONENT  ALTERNATIVES 


(1) 

(2) 

(3)  (4) 

(5) 

(6) 

(7) 

(3) 

(9) 

REPLACE 

ORG 

ORG 

ORG  ORG 

ORG 

DSU 

DSU 

DSU 

DSU 

•13PAIR 

ORU 

DSU 

GSU  Depot 

Threw 

Away 

DSU 

GSU 

Dapst 

Throw 

Away 

TABLE  3.1.  POLICY 

NUMBER  TABLE 

1  -  ORC 

2 

-  DSU 

3  -  GSU 

4 

-  DEPOT 

5  **  Throwaway 

Policy  Nuwbar  End  Itm  Rapalr  Component  Repair  Modulo  Repair 


1 

2 

1 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


18 


1 

1 

1 

1 

1 

2 

2 

2 

2 

3 

3 

3 

4 

4 

5 
2 
2 
2 
2 
3 
3 

3 

4 

4 

5 


1 

2 

3 

4 

5 
2 

3 

4 

5 

3 

4 

5 

4 

5 
5 
2 

3 

4 

5 

3 

4 

5 

4 

5 
5 
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The  following  example  shows  how  ths  fomulacor  derives  failure  node  costs 
from  cost*  for  Individual  components  and  modules. 

EXAMPLE 

Logistics  Costs  for  Component  1 

For  ORG  Removal /DSU  Repair  $100 

For  ORG  Removal /Depot  Repair  $140 


Logistics  Costs  lor  Module  1 

For  DSU  Reaoval/Throwaway  $  SO 

For  Depot  Ramoval/Throwaway  $  30 


Failure  Mode  1 

Defined  ae  Failure  of  Module  1  Causing  Failure  of  Component  1 
Accounts  for  40Z  of  Covenant  1  Failures 

Accounts  for  80%  of  Module  1  Failure  (this  implies  module  1  is  In 
sosm  other  component  besides  component  1) 


Logistics  Cost  for  Failure  Mode  1 
Policy 

OKG  Removal /DSD  Repair/Module  Throwaway 
ORG  Removal/Depot  Repair/Module  Throwaway 


Cost 

(40%) ($100)  +  (80%) ($50) 
(401) ($140)  +  (801) ($30) 


3.3  Whv  Mimed  Inteaer  Prograamlnf  la  Heeded 

The  procedure  just  outlined  develops  the  logistical  cost' Impact  of  each 
maintenance  alternative  for  each  failure  mode.  The  least  cost  set  of  alternatives 
Is  readily  identified.  If  this  set  la  implemented,  however,  test  equipment 
(and  special  repair  skills)  coats  could  be  excessive  because  no  effort  has  been 
made  to  place  all  maintenance  functions  which  use  the  same  equipment  at  as  few 
sites  as  possible.  Conversely,  if  policies  are  selected  simply  to  minimise 
test  equipment  costa,  logistical  costs  nay  become  excessive.  If  each  piece  of 
equipment  were  used  for  just  one  failure  node,  we  could  include  test  equipment 
costs  In  with  logistical  costs  In  choosing  the  policy  alternative  for  that  mode. 
When  teet  equipment  has  many  uses,  however,  we  can  no  longer  select  the  policy 
for  one  failure  mode  Independently  of  the  policies  we  choose  for  other  failure 
modes. 

A  successful  approach  to  doing  repair  level  analysis  with  shared  test 
equipment  has  been  developed  and  Implemented  by  the  US  Air  Force  building  op  work 
by  MITRE  Corporation  19,7].  This  work  does  not  Incorporate  subroutines  which  com¬ 
pute  optimum  stockage,  nor  does  it  choose  policies  subject  to  a  constraint  on 
operational  availability,  both  being  features  which  could  be  added.  Other 
restrictions  were  of  concern.  The  Air  Force  approach  appears  to  be  limited 
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to  examining  no  more  than  two  achalona  at  ona  tine.  Also,  it  doaa  not  relate 
the  maintenance  declalone  to  the  quantity  of  a  teat  equipment  required,  juat 
the  location. 

The  relatlonehip  between  maintenance  policies  and  quantity  la  not  always 
Important.  In  some  caaas,  there  la  never  a  need  to  place  more  than  one  each  of 
a  test  equipment  type  at  any  site.  In  other  caaea,  however,  maintenance  policies 
for  different  failure  modes  must  be  coordinated  not  only  to  reduce  the  number  of 
locations  for  test  equipment,  but  the  quantity  per  location;  e.g.,  it  may  be 
possible.  If  policies  are  properly  selected,  to  get  by  with  three  at  depot,  and 
one  each  at  ten  DSU's,  rather  than  one  at  depot  and  two  each  at  the  ten  DSU's. 

To  help  reduce  quantity,  0A1MBAL  may  occasionally  select  two  policy  alternatives 
for  one  failure  mode,  perhaps  depot  repair  to  handle  overflow  from  GSU  repair. 

The  Air  Force  approach  la  baaed  on  the  mathematical  technique  of  network 
theory  analyala.  Network  theory  la  a  epeclal  case  of  mixed  Integer  prograomlng 
(HIP)  In  that  any  problem  which  can  be  aolved  by  network  theory  analysis  can  be 
solved  by  HIP.  This  certainly  doea  not  work  In  reverse,  not  all  HIP  problems 
can  be  aolved  by  network  theory.  The  use  of  HIP,  as  described  In  the  next  section. 
Is  therefore  a  more  general  approach,  but  will  not  be  aa  efficient  for  problems 
where  network  theory  la  suitable. 

Interestingly,  It  la  not  oversimplifying  too  much  to  describe  most  MIP 
algorithms  as  a  synthesis  of  linear  programming  and  a  technique  called  branch 
and  bound.  Special  caaea  of  the  repair  level  analysis  problem  can  be  solved 
by  branch  and  bound  without  linear  programming  16]. 

3.4  Mathematics  of  Problem  Formulation 

The  MIP  objective  Is  to  minimise  the  sum  of  those  equipment  and  logistic 
costs  described  In  Chapter  II.  The  decision  variables  specify  where  repair  is 
to  be  done  and  the  quantity  and  placement  of  test  equipment.  The  MIP  constraints 
insure  that  all  necessary  repair  work  Is  accounted  for  and  that  the  equipment 
decisions  are  consistent  with  the  repair  decisions  In  that  the  equipment  provided 
will  handle  the  workload  placad  on  teat  equlpetent  at  each  echelon  by  the  repair 
decisions. 

Satisfying  a  system  availability  performance  goal  is  not  explicit  in  the 
objective  functions  or  constraints,  but  does  Impact  on  the  policies  chosen.  How 
this  Is  accomplished  Is  explained  In  Sections  3.5  and  3.6. 

First,  let's  discuss  the  constraint  and  objective  function  equations  of  the 
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HIP  In  nor*  detail.  Notation  and  tha  entire  formulation  are  suamarlsad  in  Section 
3.3.  Wa  will  be  referring  to  the  equations  of  Section  3.3,  one  equation  at  a  time, 
beginning  with  the  conatralnta  and  then  examining  the  objective  function. 

Every  t late  a  failure  node  Incident  "Ij"  occurs  (nodule  j  In  component  1 
fella  ao  the  end  It  me  Is  repaired),  maintenance  actions  are  generated  and  some 
maintenance  policy  nuat  react  to  than.  For  example,  the  end  item  mould  be  re¬ 
paired  at  echv ion  ?  by  replacing  component  i;  the  component  would  be  repaired 
at  echelon  R  by  replacing  module  j;  and  the  module  would  be  repaired  at  level 
r.  A  decision  to  be  made  by  the  MIP  la  what  percentage  of  ouch  failure 
incidents  should  be  handled  by  policy  (P,Rtt),  and  the  constraint  (for  each 
failure  node)  la  that  the  percentages  sum  to  1  over  all  feasible  policies. 

Kchelona  are  numbered  ae  follows: 


Echelon 

ORG 

DSU 

GSU 

Depot 

Throwaway 

Number 

1 

2 

3 

4 

5 

Currently,  P  must  be  either  1  or  2  (end  item  repair  at  ORG  or  DSU);  R  must  equal 
or  exceed  P  (component  cannot  be  repaired  at  lower  echelon  than  end  item);  r 
must  equal  or  exceed  R  (repair  module  at  or  above  where  it  is  replaced) .  In 
Section  3.S,  the  constraint  on  percentages  Is  denoted  the  "Accountability" 
constraint,  all  work  la  accounted  for. 

The  next  set  of  constraints  relate  to  workload  on  test  equipments.  Workload 
la  input  to  OATMEAL  aa  the  variables  N^(e) ,  n^  (e) ,  N^(e),  depending  on  the 
function.  N^(a)  defines  workload  for  equipment  e  attributable  to  repair  of  the 
end  Item  using  component  1,  so  that  any  failure  mode  Involving  component  i  places 
a  workload  of  N^e)  on  teat  equipment  e.  n^(e)  defines  workload  attributable  to 
the  repelr  of  module  j,  while  (e)  defines  workload  attributable  to  repair  of 
component  1  when  module  j  falls.  Workload  is  defined  ss  the  fraction  of  equipment 
availability  required  per  end  item  per  year:  if  equipment  e  is  up  2000  nours 
a  year,  if  repair  of  the  end  item  when  component  1  fails  requires  2  hours  use,  and 
if  the  component  1  failure  rate  is  0.5  times  per  year  per  end  item,  then  Ni(e) 
equals  (2.0  x  0.5)/2000  or  .0005. 

There  is  a  workload  constraint  for  each  equipment,  e,  for  each  echelon  k. 

If  Ni(e)  is  non-zero,  any  failure  mode  involving  component  i  with  end  item 
replacement  at  echelon  k  contributes  to  the  workload;  put  differently,  if  Y 
(k,R,r)  Is  non-zero  for  some  values  of  R  and  r,  and  N^(e)  is  non-zero,  there  is 
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a  contribution  of  [Y^ (k,R,r)]  x  iN^CeJJ  x  Iequipment  density  per  echelon  k  site] 
to  the  workload  constraint  for  equipment  e  at  echelon  k.  Similarly,  if  nj(e)  and 
Y^j(P,R,k)  are  non-zero  for  any  values  of  P  and  R  there  is  a  contribution  to 
workload  at  echelon  k  as  there  is  if  N^ (e)  and  Y. j  (P,k,r)  are  non-zero. 

The  sum  of  all  3  types  of  requirements  attributable  to  N^e),  N^e),  n^  (e) 

must  be  less  than  or  equal  to  t  .  ,  the  number  of  test  equipment  e  at  echelon  k 

e(  k 

in  order  to  Insure  (in  a  steady  atate  sense)  that  the  repair  functions  can  be 
performed.  These  t  .  variables  are  pseudo-decision  variables  in  that  they  are 

6|  K 

almost  entirely  dependent  upon  the  Y^  decisions;  in  the  cases  where  the  number 
of  test  equipment  e  at  an  echelon  k  must  be  an  Integer  value  (TE  is  peculiar  to 
this  end  item)  there  is  interdependency  in  that  the  MIP  procedure  may  modify 
the  Y. .  decisions  to  minimize  integers  t  .  . 

The  objective  function  is  to  minimize  costs  associated  with  the  Y^.,  k 

decisions.  The  te  ^  values  are  multiplied  by  the  equipment  cost  C@  anu  the 

number  of  repair  sites  at  echelon  k.  Associated  with  the  Y^CP^r)  decisions 
are  the  logistics  costs  for  the  failure  mode  "ij"  and  policy  (P,R,r).  Section  3.1 
discussed  computation  of  logistics  coats  by  failure  mode  in  detail.  Summarized 
mathematically,  the  component  and  module  logistics  costs,  C^PjR)  and  Mj(F,r) 
respectively,  are  prorated  by  fractions  FC^Cj )  and  (i)  which  denote, 
respectively,  fraction  of  component  i  removals  accounted  for  by  failure  mode 
"ij"  and  fraction  of  module  J  removals  accounted  for  by  failure  mode  "ij." 

Size  of  MIP.  There  are  twenty-five  possible  combinations  of  values  for 
(P,R,r) :  (1,1,1),  (1,1,2),  (1,1,3),  (1,1,4),  (1,1,5),  (1,2,2),  (1,2,3),  (1,2,4), 

(1,2,5),  (1,3,3),  (1,3,4),  (1,3,5),  (1,4,4),  (1,4,5),  (1,5,-),  (2,2,2) . (2,5,-). 

Let 

NAPP  -  number  of  failure  modes,  sometimes  referred  to  as  applications. 

NEQECH  -  number  of  equipment/echelon  combinations,  so  that  if  there  are 
10  equipment  types,  each  of  which  can  be  placed  at  any  of 
3  echelons,  NEQECH  ■  30. 

Then: 

Numer  of  Accountability  Contraints  1  NAPP 

Number  of  Workload  Constraints  *  NEQECH 

Number  of  Continuous  Variables,  Y^(P,R,r)  ■  25  x  NAP? 

Number  of  possibly  integer  variables  t  ,  ■  NEQECH. 
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3,5  Mathematical  Formalism  of  ORLA  Model 


Inputs  to  ORLA 
C^(P,R)  - 

Mj (R,r)  - 


DENS, 

k 


nj <e)  ' 


(e)  - 


Nt(e) 


F^O)  - 

FMj (i)  - 


total  repair  and  logistics  costs  associated  with  component 
i  vhen  used  at  level  P  (to  fix  the  end  item)  and  repaired 
at  level  R. 

total  repair  and  logistics  costs  associated  with  module  j 
when  used  at  level  R  (to  fix  a  component)  and  repaired 
at  level  r. 

density  of  equipment  supported  by  an  echelon  k  supply/repair 
facility. 

fraction  of  1  year's  working  hours  of  equipment  e  required 
to  repair  module  j  for  all  failures  of  module.  (Per  end  item) . 
fraction  of  1  year's  working  hours  of  equipment  e  required 
to  repair  component  i  when  module  j  falls  for  all  failures 
of  that  mode.  (Per  end  item). 

fraction  of  1  year's  working  hours  of  equipment  e  required 
for  repair  of  the  end  item  for  all  failures  of  component  i. 

(Per  end  item). 

number  of  echelon  k  supply/repair  facilities. 

cost  of  equipment  or  mos  type  e . 

percent  of  component  1  failures  due  to  module  j . 

percent  of  module  j  failures  which  occur  in  its  application 
to  component  1. 


Decision  Variables 

t  ,  -  number  of  equipment  e  at  echelon  k  (defined  as  Integer  for 

e  j  K 

peculiar  test  equipment /repair  skills). 

(P,R. r)  -  A  failure  mode  is  designated  by  the  component  (i)  and 

module  (J)  involved.  This  variable  gives  the  percent  of 
failures  for  that  mode  for  which  the  policy  is  to  repair 
the  end  item  at  echelon  P,  repair  the  component  at  echelon 
R  and  repair  the  module  at  echelon  r.  If  P  is  5,  this 
means  component  is  thrown  out.  If  R  is  5  it  means  the 
module  is  thrown  out. 


Minimize 


+ 


2 


H 


ii 


Subject  to: 


E  C 


kIUk  ‘e.k 


4*  5  5  r 

E  1  E  r  E  Y  <P,R,r)  *  |[FC.(j)L(P,R)  + 

i  j  P-1  R-P  r-R  L 


[FM  (DKMj  (R.r)) 


} 


4  5 

E  Z 

P-1  R-P 


DENS.  E 
k  i 


+  DENS,  E 
k  1 


+  DENS,  E 
k  1 


0 


5 

E  Y  <P,R,r)  -  1 
r-R 


"Accountability" 


"Availabilicy" 


All  Variables  1  0 

t  ,  integer  if  the  equipment  is  peculiar 

e,  k 


is  < 


1. 


3.6  MIP  Approximation 

Aa  mentioned  In  Section  3.1,  the  method  cf  computing  logistics  costs  by 
failure  mode  using  proretlon  la  only  approximately  correct.  For  example,  suppose 
a  component  has  two  failure  modes,  and  that  the  MIP  selects  a  policy  of  DSU  re¬ 
pair  when  failure  mode  1  occurs  and  depot  repair  when  failure  mode  2  occurs. 

Let  the  facts  be: 

Logistics  Costs  for  component 

for  1001  DSU  repair  -  $100 
for  1001  depot  repair  -  $140. 

Percent  of  component  failures  due  to  each  failure  mode 
failures  due  to  mode  1  -  60% 
failures  due  to  mode  2  -  401 

The  total  logistics  costs  for  the  component  as  computed  In  the  MIP  objective 
function  would  be  601  ($100)  +  401  ($140). 

A  more  accurate  assessment  of  costs  would  be  obtained  by  running  SESAME 
with  a  Maintenance  Task  Distribution  showing  601  DSU  repair  and  401  depot 
repair.  In  fact  this  would  be  done  In  the  EVALUATOR.  This  cannot  be  done  in 
the  MIP  FORMULATOR  which  develops  the  objective  function,  because  it  is  only 
after  the  MIP  is  solved  that  we  know  what  will  be  done  for  each  failure  mode 
and  therefore  what  the  Maintenance  Task  Distribution  should  be. 


3.7  Relation  of  Supply  Lagranglan  to  Maintenance  Policy 


For  a  given  maintenance  policy  vector  (a  designation  for  every  component  and 
module  in  the  system  of  where  to  replace  and  repair) ,  a  definite  relation  exists 
between  the  SRSAME  Lagranglan  value  and  the  availability  achieved  (higher  value, 
higher  availability);  and  this  relation  is  developed  thru  the  SESAME  multi-echelon 
stockage  decisions.  However,  beforv  the  policy  vector  Is  fixed  (by  the  OATMEAL 
optimization  procedure)  the  relation  between  a  chosen  Lagranglan  value  and  an 
achievable  system  availability  Is  not  necessarily  a  one  to  one  mapping. 

For  example  e  maintenance  policy  vector  P(l)  la  initially  chosen;  a 
Lagranglan  value  *T(1)  Is  found  that  compels  enough  stockage  so  that  the  target 
availability  A^  Is  achieved.  For  that  1^(1),  the  MIP  finds  a  maintenance  policy 
P(2)  that  minimizes  stockage  backorder  end  maintenance  costs  but  achieves  an 
availability  Ap^*  In  order  to  now  meet  A^  with  policy  vector  P(2)  the  Lagranglan 
may  have  to  be  raised  or  lowered  to  X^(2).  However*  now  the  policy  P(2)  may  no 
longer  be  optimal  for  X^(2).  More  about  this  is  detailed  in  Appendix  A.l. 

To  avoid  this  potential  for  looping  a  method  for  Intelligently  choosing 
an  initial  Lagranglan  Is  needed*  but  let's  first  review  the  whole  Lagranglan 
concept . 

3.8  initial  Lagranglan  Selection  and  Treatment  of  Operational  Availability 

Constraint 

Using  the  formula  for  operational  availability  (OA)  in  Section  2.6,  a  con¬ 
straint  on  OA  may  be  stated  as: 


MCTBF _ 

MCTBF  +  MTR  t-  MTT  +  MLDT  - 


TARGET 


By  algebraic  manipulation,  this  Is  equivalent  to: 


MCTBF  -  (MCTBF) (TARGET) 
TARGET  - 


MTR  +  MTT  +  MLDT 


The  repair  level  analysis  objective  may  then  be  stated  as: 


Minimize  Cost 


Subject  To:  MTR  +  MTT  +  MLDT  —  (ftCTBF)  (TARGET} 

TARwlST 


Under  the  Generalized  Lagranglan  method,  this  problem  is  transformed  to 

an  Unconstrained  optimization  with  a  Lagranglan  parameter  "X": 

26 


Minimize  Coat  +  (A)  (MTR  +  MTT  +  MLDT) 


A  solution  to  this  problaa  will  have  associated  with  It  both  a  cost,  and 
an  achieved  OA.  Everett’s  Theorem  {2]  guarantees  that  no  other  solution  with 
equal  or  higher  operational  availability  can  cost  less. 

Thus,  one  approach  to  acconodatlng  the  OA  constraint  is  depicted  In 
Figure  3.1:  we  choose  a  Lagrangian,  solve  the  transformed  problem  and  adjust 
the  Lagrangian  until  we  find  a  solution  with  OA  close  to  our  target. 

Recall  that  one  set  of  variables  in  the  MIP  objective  function  are  the 
C^P.R).  To  solve  the  Lagrangian  fora  of  the  problem,  the  end  item  delay  associated 
with  replacing  component  i  at  level  P  a-*d  fixing  it  at  level  R  is  multiplied  by 
A  and  added  to  (^(P.R).  Delay  will  be  caused  either  by  component  backorders  - 
this  value  and  its  Impact  on  MLDT  is  computed  by  the  SESAME  subroutines  used  in 
getting  the  C^(P,R)  -  or  because  P  is  above  ORG,  adding  to  MTT  (cf  Section  2.6). 

The  difficulty  with  the  approach  as  outlined  is  that  it  can  be  time  consuming 
since  the  whole  problem,  using  MIP,  must  be  resolved  each  time  a  new  value  of  A 
is  tried. 

For  many  applications  the  user  is  willing  to  accept  some  degree  of  non- 
optimality  in  order  to  have  a  tool  which  is  easy  to  use.  We  will  therefore  consider 
how  to  avoid  looping  and  still  obtain  a  good  answer. 

What  we  would  like  to  do  la  find  a  good  initial  value  for  "A".  Then  once 
the  MIP  is  run,  any  discrepancy  between  achieved  and  target  OA  is  addressed  by 
modifying  stockage  policies,  while  retaining  the  maintenance  policies  found  in 
the  MIP.  The  Evaluator  is  programmed  to  do  this;  given  the  set  of  maintenance 
policies  found  by  the  MIP,  it  will  always  compute  stockage  based  on  the  target  OA; 
it  does  not  simply  reproduce  the  stockage  quantities  Incorporated  in  the  C^(R,P) 
during  the  MIP  formulation. 

We  would  like  to  determine  as  accurately  as  possible,  before  running  the  MIP, 
the  relationship  between  the  Lagrangian  value  used  and  the  resulting  OA  property 
of  the  MIP  solution.  This  relationship  depends  on  where  each  component  Is  replaced 
and  repaired  in  the  MIP  solution. 

Most  critical  is  where  the  component  is  replaced.  This  completely  determines 
the  contribution  of  that  component  to  MTT.  (Section  2.6).  The  repair  echelon 
helps  determine  the  contribution  of  that  component  to  MLDT,  but  it  is  not  critical 
to  the  value  of  MLDT  which  will  emerge.  If  a  component  is  repaired  at  depot  rather 
than  DSU,  more  stockage  may  be  required  to  cover  a  longer  pipeline,  but  the 
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contribution  to  MLDT  need  not  change  much.  In  fact  the  fill  rate  for  a  component 
at  user  level  in  the  optimal  solution  depends  "essentially"  only  on  the  Lagrangian 
value,  not  on  where  the  component  is  repaired.  This  is  discussed  in  [3]  where 
sufficient  background  Is  given  to  explain  the  Impact  of  "essentially." 

These  observations  were  Incorporated  into  the  method  for  choosing  the 
initial  "1".  A  search  routine  for  X  is  built  into  the  MIP  formulation;  it  is 
comparable  to  the  search  routine  incorporated  into  the  full  SESAME  supply  model 
which  is  also  designed  to  find  the  X  which  will  give  the  target  OA.  The 
difference  is  that  in  SESAME  ve  know  in  advance  the  replace  and  repair  echelons 
for  each  component.  In  the  MIP  formulator  we  guess  that  the  repair  echelon  will 
be  the  GSU  since  this  Is  not  critical,  and  determine  the  replace  echelon  as  part 
of  the  search. 

The  whole  process  is  depicted  in  Figure  3.2.  For  each  trial  x*  stockage  is 
computed  using  the  appropriate  subroutines  from  SESAME,  and  the  resulting  contri¬ 
bution  of  that  component  to  MLDT  and  MTT  is  determined  by  these  subroutines. 

Once  all  items  have  been  processed  the  OA  can  be  computed  and  x  adjusted  as 
necessary.  In  the  MIP  formulator  each  item  is  processed  twice,  once  assuming 
the  component  is  replaced  at  O&G,  and  once  assuming  it  is  replaced  at  DSU. 

How  io  we  know  which  set  of  answers  to  use,  those  based  on  the  ORG  replace¬ 
ment  assumption  or  those  basec*.  on  the  the  DSU  assumption?  For  each  component 
this  is  based  on  the  replacement  assumption  which  leads  to  least  cost  for  that 
component.  Coat  is  based  both  on  logistical  cost,  and  on  a  proportional  share 
of  potential  test  equipment  coot.  For  example,  in  assessing  cost  given  ORG 
replacement,  we  assume  all  teat  equipment  needed  for  end  item  repair  is  at  ORG  and 
allocate  cost  to  the  component  baaed  on  the  test  equipment  throughout  it  requires  as 
a  fraction  of  total  requirements  for  that  test  equipment  for  end  item  reoair. 

To  summarize,  initial  x  selection  is  based  on  a  search  routine  in  the  MIP 
formulator,  akin  to  the  search  for  "x"  in  the  SESAME  supply  model.  In  this  search 
we  try  to  predict  the  replacement  level  for  each  component  which  will  be  chosen 
by  the  MIP.  We  will  not  always  be  right  because  our  way  of  handling  test  equip¬ 
ment  costs  in  the  search  is  not  comparable  to  the  way  the  MIP  does  it. 

Experience  with  test  cases  on  the  x~*eerch  procedure  is  discussed  in 
Appendix  A. 2.  Appendix  A. 3  discusses  computer  costs,  and  how  to  proceed  In 
those  cases  where  the  user  wishes  to  devote  the  effort  to  getting  the  ideal  value 
for  x  as  depicted  in  Figure  3.1. 


FIGURE  3.1 


APPENDIX  A. 


INPUTS  TO  OATMEAL 

A.  1  Coat  Parameters  and  Operational  Availability  Data 

The  axpraaalona  on  tha  right  art  as  defined  In  Chapter  2.  Numbers  In 
parentheala  indicate  tha  variable  ia  a  vactor  with  number  of  elements  as  specified. 


COSNSN 

COSBIN 

COSREQ 

C0STRA(3) 

PVPAC 

CHFAC 

AVTAR 

MTR 

MCTBF 

CTDFL 


Item  Introduction  Coat  +  PVPAC  x  Item  Maintenance  Cost 

PVPAC  x  Bln  Coat  Kate 

PVPAC  x  Requisition  Coat  Rata 

* 

Transportation  Rates  :  1  la  DSU  to  ORG,  2  is  GSU  to  DSU, 

3  la  Depot  to  GSU 

PVPAC 

(1.  +  HCFAC  x  PVPAC) 

Operational  Availability  Target;  e.g.  0.95 

Mean  Tima  to  Repair ,  in  hours 

Mean  Calendar  Tima  Between  Failures,  in  days 

"Contact  Team  Delay."  Tima  lost  when  end  item  repair  is  not 
accomplished  with  ORG  level  resources,  in  days. 


A.2  Supply  System  Data 

These  variables  are  all  as  defined  in  the  SESAME  User's  Guide  [10],  end  item 
card:  CIMNTS(3),  OST(4),  OPSL(A),  SSC,  KSC.  World  wide  density  is  input,  and  all 
OUPS  are  calculated  from  this  and  number  of  claimants  (OUPS  »  density  i  claimants). 


A. 3  Maintenance  System  Data 

NEQ  Number  of  Test  Equipments 

NLRU  Number  of  Components,  also  labelled  "LRU's" 

NSRU  Numer  of  Modules,  also  labelled  "SRU's" 

NAPP  Number  of  Failure  Modes,  also  labelled  "applications" 

ERPATE  Error  or  false  removal  rate.  If  ERRATE  is  10%,  failure 
rates  are  multiplied  by  1.1  to  get  removal  rates. 

TERAT(A)  Ratio  of  work  week  hours  for  each  echelon  to  hours  at 

echelon  1.  Thus  TERAT  (1)  is  always  1.  This  is  used  in 
computing  requirements  for  test  equipment. 


A.  4  Policy  Constraints 

A  maintenance  policy  as  it  pertains  to  a  given  failure  mode  can  be  expressed 
as  a  triplet,  (i,j,k),  where  i,j  and  k  specify  where  respectively  the  end  item, 
a  _  > 

When  Supply  Structure  Option  Code  is  D  or  N,  and  there  is  direct  ordering  from 
Depot  to  DSU,  OATMEAL  will  use  max  [C0STRA(2),  C0STRA(3) ] .  C0STRA(2)  will  be 
applied  to  retrograde  from  DSU  to  GSU. 


component ,  and  module  art  repaired.  Since  there  are  only  four  echelons,  a  5 
denotea  throwaway.  Table  4.1  aaaoclataa  a  number  with  each  policy  triplet. 

The  uaer  inputs  the  array  IPOL.  Only,  if  IPOL (a),  a  •  1  to  25  is  1, 
le  policy  a  a  candidate;  thus,  if  IPOL (3)  le  0,  the  policy  triplet  (1,1,3) 
la  not  allowed. 


A.  5  Teat  Equipment  Pets 


The  following  are 
EQCS?(4) 

IEQPEC 


IEQPLA 

EQSTK(IOOO) 


input  for  each  piece  of  teat  equipment: 

[Purchaae  Coat  +  Maintenance  Coat  a  PVFAC)  as  defined  in 
Section  2.4. 

The  equipment  la  coanon  at  all  echelons  above  IEQPEC: 

If  IEQPEC  is  0,  it  la  totally  common,  if  IEQPEC  la  4, 
it  la  totally  peculiar. 

The  equipment  nay  only  be  placed  at  echelon  IEQPLA  or 
above. 

Each  entry  conslata  of  a  decimal  of  the  form  n.x.  "n" 
is  an  equipment  number,  "x"  is  the  fraction  of  that 
equipment' a  annual  throughput  used  per  maintenance  action. 
Entries  are  aasociated  with  specific  components,  modules 
and  applications.  If  an  entry  with  value  4.01  is  associated 
with  a  component  it  means  that  when  that  component  is  used 
to  fix  a  system,  test  equipment  4  is  needed,  and  100  such 
maintenance  actions  per  year  could  be  handled  by  equipment 
4  if  it  were  devoted  exclusively  to  that  type  action.  If 
an  entry  la  associated  with  a  module,  it  refers  to  mainte¬ 
nance  of  the  module.  If  an  entry  is  associated  with  an 
application,  it  refers  to  repair  of  the  component  using 
the  module. 


A. 6  Component  and  Module  Data 

UP+  Onit  Fries  x  CHFAC 

WGT  Weight  in  Pounds 

IESS  Essentiality  Code 

WASU  Percent  of  removed  items  which  will  washout  as  unrepair¬ 

able  regardless  of  policy. 

TAT(4)  Shop  turn  around  time  in  days  to  fix  module  or  component 

by  echelon. 

32 


■v  "k 


X  \ 


INDSTK  Lut  entry  In  EQ8TK  (defined  above)  aaeoclated  with  the 

item. 

N  STACK  Number  of  entriee  in  FQSTK  aaeociated  with  the  item. 

PARTSR  Number  of  nee  parts  per  module  (0  for  conponent) 

PARTS P  Average  price  of  parte 

DOC  "One  Tine  Repair  Costs,"  euch  aa  manual  preparation 

(documentation)  coete.  (See  Section  2.4.2) 

REPC(I)  The  "Repair  Coat^"  (See  Section  2.4.2) 

NNSN  If  0,  item  is  new  to  syatam;  fractional  value  can  be  used 

when  RN  >  1. 

ID  Alphanumeric  ID 

RN  Repetition  Nuaber  (Used  when  an  "average"  component  or 

module  represents  a  nuaber,  RN,  of  components  or 
modules  not  entered  individually) . 

A.  7  Failure  Mode  Application  Date 

Identification  nuaber  of  component  to  which  failure  mode 
pertains.  A  component  or  module  number  is  determined 
by  the  order  in  which  it  is  read  in;  e.g.  the  5th 
component . 

Identification  nuaber  of  aodule. 

* 

Number  of  failures  per  end  item  per  year. 

Shop  turn  around  time  to  fix  component  if  module  fails. 

Reference  entries  in  EQSTK. 

A.  8  Inputs  to  SESAME  SUBROUTINES  of  OATMEAL 
A. 8.1  CURPARs 

For  components,  the  Lagranglan  is  used.  (See  Sections  3.5  and  3.6). 
For  modules,  the  price  of  the  component  on  which  it  appears  is  used,  or  a  weighted 
average  if  the  module  appears  on  more  than  one  component.  For  parts,  the  price 
of  the  module  is  used.  This  is  consistent  with  how  part  and  module  backorders 
are  costed  out  (Section  2.4.9). 

A.8.2  WHOFIL  and  CONDEL 

Wholesale  stockage  is  computed  for  each  item,  and  the  expected  W110FII. 
and  CONDEL  which  will  result  from  this  stockage  are  used  for  that  item  in  com¬ 
puting  retail  stockage;  thus  WHOFIL  and  CGNDEL  vary  by  item  and  are  consistent 
with  the  dollar  cost  of  wholesale  stockage. 

*Whan  FAIL  pertains  to  average  component  or  module  (RN  greater  than  1) ,  t All 
would  be  multiplied  by  RN  to  get  total  failures. 


IDL 


IDS 
FAIL 
TAT  (4) 
NSTACK 
INDSTK 
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Wholesale  atockage  la  baaad  on  a  raorder  quantity  of  1.  The  reorder 
point  la  baaad  on  a  target  atock  availability  computed  aa: 


“  (.60  ♦ 


COTTAR 
COTTAR  +  UP; 


Thla  la  an  average  of  .60  and  the  availability  target  appropriate  at  uaer  level. 
It  la  conalatent  with  the  findings  on  optima  upper  echelon  availability  targeta 
documented  in  [3]. 

A. 8. 3  TATa  and  Rotrotrado 

The  uaer  Inputa  ahop  repair  time.  Turn  around  tlmea  Input  to  SESAME 
include  retrograde  tlmea.  It  la  aaaumad  retrograde  tlmea  between  the  echelons 
equal  order  and  ahlp  tlmea  between  thoae  echelona.  It  la  alao  aaaumad  the  lower 
the  echelon  at  which  a  repairable  la  removed,  the  lower  the  echelon  at  which 
It  la  repaired. 


Example: 

ORG 

DSU 

GSU 

WASHOUT 

RTD 

60X 

AOX 

MID 

AOX 

30X 

30X 

Assumption: 

All  ORG  removals  are  fixed  at  DSU  (AOX  of  60Z)  and  GSU  (201  of  60X). 

DSU  removals  are  fixed  at  GSU  (10X  of  40X)  and  are  washed  out 
(30X  of  AOX). 

Finally,  In  a  non-vertical  system,  where  the  GSU  Is  repairing  items 
for  the  DSU  stock  account,  turn  around  time  will  include  the  time  to  get  the  good 
Item  from  GSU  to  DSU. 

The  sane  assumptions,  in  fact  the  same  computer  subroutine,  used  to 
get  retrograde  times  are  alao  used  to  get  retrograde  coats.  Instead  of  time  being 
assessed  for  each  retrograded  item  (and  added  to  repair  time) ,  a  coat  is  assessed 
based  on  the  retrograde  coat  rate. 

A. 8. A  Replacement  Task  Diatribution  and  Maintenance  Task  Distributions 

(RTD.  MTpy  “ 

An  RTD  Is  a  percentage  breakout  of  the  replacements  of  a  component 
or  module  by  echelon;  similarly  an  MTD  la  a  percentage  breakout  of  repair  of  a 
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component  or  module  by  echelon,  Including  the  throwavay/vashout  percent 
referenced  by  echelon  "5".  It  le  part  of  the  EVALUATOR* a  job  to  build  the 
RTD^,  KTD^  of  a  component  i  and  the  RTD, ,  HTD^  of  a  nodule  j  from  the  policy 
variable  valuee  of  the  failure  nodea  that  pertain  to  component  J  or  module  j. 

Remember  that, 

*ij(*»R»r)  ■  percent  of  failurea  of  nodule  J  in  component  1  for  which 
policy  le  to  replace  i  at  level  ?,  repair  1  by  replacing  j  at  level  F  and  repair 
j  at  level  r. 

Ualng  theae  value*  and  the  proportion  of  failures  of  component  1 
due  to  module  J  and  the  proportion  of  module  j  failurea  that  occur  in  component  i, 
one  can  build  up  the  percentage*,  by  echelon,  of  replace  and  repair  associated 
with  P,  R  and  r  type  decieiona  for  component  i  and  module  j.  Note  that  in  this 
procedure  one  must  adjuat  the  module  RTD'e  and  MTD'a  for  instances  when  the  next 
higher  assembly  component  wee  thrown  out. 

The  eggregetlon  or  buildup  of  teek  distributions  is  done  by  summations 
over  the  proper  indices.  For  example,  the  RTD  for  component  i  at  echelon  is 
found  by  e  suomation  process  over  j,R,r  of  instances  involving  index  i,  and  I*  -  k. 

A.  9  Promotion  of  Modulee 

If  the  component  le  repaired  at  the  same  echelon  as  the  end  item  und  if 
there  is  1001  repair  (no  washout),  it  is  assumed  the  component  will  net  be 
stocked;  Instead  it  will  be  repaired  eo  that  the  same  component  remains  in  the 
end  item.  This  situation  has  the  following  ramifications: 

a.  Repair  time  for  the  component  directly  degrades  operational 
availability. 

b.  Backorders  of  a  module  used  to  repair  the  component  lengthen  com¬ 
ponent  repair  time  and  so  degrade  operational  availability.  The  module  is 
"promoted,"  and  must  be  treated  as  an  LRU. 

In  the  HIP  Formulator,  when  an  estimate  of  CURPAR  is  derived,  the  possibility 
>f  module  promotion  is  not  considered.  Otherwise,  all  ramifications  of  iroduK- 
promotion  are  Implemented. 


APPENDIX  B 


PROPERTIES  OF  LAGRANGIAN  AND  MAINTENANCE  POLICIES 


Relation  of  Sut 


Lag rang ian  to  Maintenance  Polic\ 


Figure  4  depicts  a  family  of  availability  vs  cost  curves .  The  horizontal 
axis  consists  of  maintenance,  transportation,  supply  costs  -  i.e.,  all  those 
costs  in  the  ORLA  model,  excluding  backorder  penalty  costs,  that  would  be  affected 
by  a  particular  maintenance  posture  -  and  the  vertical  axis  nipresenus  the  avail¬ 
ability  achieved  as  the  Lagrangian  penalty  cost  (SESAME  curve  parameter  CURPAR) 
is  varied  and  stockage  costs  increased.  Each  single  curve  represents  one  main¬ 
tenance  policy  vector  -  a  designation  for  every  LRU  and  module  in  the  system  of 
where  to  replace  and  repair  components.  The  number  of  combinations  of  designa¬ 
tions  could  be  enormous  and  the  figure  depicts  this  with  a  partial  representation 
of  closely  spaced  curve  spectra.  As  one  moves  rightward,  the  curves  represent 
increased  use  of  fix  forward  (more  parts  fixed  at  lower  echelon);  repairing  at 
ORG  and  DSU  incurs  relatively  higher  costs  even  for  low  availability  (more  TE, 
personnel,  stockage  quantities)  but  the  potential  for  quick  repair  response 
allows  higher  availabilities  to  be  attainable.  Conversely,  for  curves  beginning 
at  the  lefthand  portion  (more  and  more  parts  maintained  and  stocked  at  higher 
echelons  -  GSU,  DEPOT)  fixed  costs  are  generally  lower  but  availability  can  be 
increased  (by  stocking)  only  so  much  due  to  inherent  transport  delays  of  end  items 
or  unserviceable  assemblies  rearward. 

The  spectral  properties  of  the  curves  would  not  be  as  pure  as  shown 
(for  clarity);  there  would  be  crossing  of  curves  due  to  transportation  and 
maintenance  cost  tradeoffs.  The  figure  can  be  considered  a  response  surface 
on  which  one  moves  to  find  the  best  solution  -  achieving  a  target  system  avail¬ 
ability  while  minimizing  the  total  of  all  costs.  The  tools  for  conducting  the 
search  are  a  Lagrangian  (tied  to  a  multi-echelon  stockage  optimizer)  to  vary 
stockage  quantities  and  availabilities,  and  a  MIP  to  find  the  maintenance  policy 
vector  that  minimizes  costs,  including  the  Lagrangian  tied  backorder  cost.  In 
other  words,  so  that  availability  is  considered  in  selecting  policies,  user  level 
component  backorders,  and  delays  caused  by  repair  of  the  end  item  above  user  level, 
are  costed  out  using  the  Lagrangian  and  included  in  logistics  costs. 

In  the  discussion  that  follows,  for  consistency  in  the  schematic  exemplifica¬ 
tion,  its’s  indicated  that  the  MIP  chooses  the  leftmost  curve  point  for  a  given 
Lagrangian  value  -  i.e.  that  which  has  minimum  cost,  excluding  backorder  costs. 
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The  same  argument  below  would  follow  If  other  points  were  Indicated  on  the  figure 
as  the  MIP  choices. 

Suppose  a  target  availability  is  chosen  to  be  A^,  as  shown  in  tuo  diagram. 

To  find  a  Lagrangian  one  can  choose  an  Initial  policy  vector,  say  curve  a  iep re¬ 
senting  the  most  fix  forward  concept;  the  choice  for  the  Lagrangian  would  be  o, 
since  this  value  on  curve  a  yields  an  availability  near  A^,.  The  MIP  would  then 
select  curve  b  for  the  final  policy  vector  (since  its  o  is  left  most)  but  evalua¬ 
tion  of  the  performance  for  that  policy  of  rear  repair  would  show  that  o  or  any 
value  higher  (a*x,q)  could  never  reach  the  target  availability.  If  one  chose 
the  next  most  right  curve,  c,  the  Lagrangian  would  be  A,  and  the  MIP  would  home 
in  on  curve  d  for  that&  ;  then  to  reach  target  the  penalty  parameter  would  have 
to  be  raised  during  an  evaluation  to  a  value  between  x  and  a,  but  at  that  point 
the  policy  is  no  longer  optimal!  Similary  guessing  policy  curve  d  and  thereby 
□,  the  MIP  would  settle  on  curve  c,  and  moving  down  that  curve  to  around  pointy  , 
indicates  also  that  that  policy  is  no  longer  the  optimal  one  in  that  region  near 
the  target. 

It  became  clear  that  we  needed  a  procedure  for  choosing  an  initial  policy 
or  Lagrangian  value  such  that  the  policy  -  Lagrangian  value  was  in  the  neighbor¬ 
hood  of  the  optimum  on  our  response  surface.  In  the  figure’s  example  a  choice  of 
curve  d  or  curve  e,  which  would  lead  to  a  chosen  value  near  x,  would  then  lead  to 
a  MIP  policy  curve  that  does  not  necessitate  much  if  any  change  in  the  Lagrangian 
to  achieve  availability  -  and  hence  no  need  to  perturb  the  final  policy  vector. 

Properties  of  the  Optimal  Solution 

Although  Figure  4  is  misleading  in  terms  of  generalizing  from  a  pedagogical 
example,  there  may  be,  as  indicated  by  x,  ,  and  o  points  in  the  region  of  target 
levels,  several  combinations  of  Lagrangian  and  policy  vector  that  can  lead  to 
the  availability  at  near  minimum  cost.  And  in  terms  of  the  cost  -  Lagrangian 
policy  surface,  it  is  probably  flat  around  the  optimum,  i.e.  sub-optimal  policies 
might  be  quite  satisfactory  in  terms  of  relative  cost  savings.  So  even  if  our 
policy-vector  found  by  the  procedure  3.6  and  Appendix  C  is  not  purely  robust,  it  is 
satisficing.  Also  we  may  allow  the  decision  maker  to  make  adjustments  to  the 
policy  on  an  individual  component  basis;  the  evaluator  will  have  a  mode  to  assess 
perturbations  to  our  final  policy. 
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APPENDIX  C 


TEST  EXPERIENCE  WITH  THE 


SEARCH 


Three  real  world  end  items  were  investigated,  an  expensive  tactical  radio 
with  counter- jamming  features,  a  simple  radar  and  a  more  complex  radar.  Table 
C.l  characterizes  these  end  items  in  terms  of  number  of  components,  modules, 
failure  modes,  and  test  equipments. 

We  were  interested  in  how  well  the  heuristic  X  selection  described  in 
Section  3,6  would  work.  Our  criteria  was  the  cost  of  the  solution  found  as 
compared  to  a  solution  found  by  looping,  also  described  in  Section  3.6.  While 
doing  this  work  we  encountered  an  unexpected  problem  with  the  MIP  procedure  it¬ 
self. 

A  MIP  procedure  has  two  objectives  -  finding  an  optimum  solution,  and 
guaranteeing  that  the  solution  found  is  truly  optimum..  To  save  running  time, 
any  solution  which  is  guaranteed  to  be  within  some  percent  or  dollar  value  of 
the  best  possible  solution  is  usually  accepted  as  optimum.  We  used  a  criteria 
of  0.5%,  but  also  set  a  maximum  on  the  amount  of  computer  resources  a  KIP  run 
was  allowed  to  use,  about  $45. 

In  the  second  iteration  of  the  loop  procedure  for  the  tactical  radio,  the 
MIP  could  not  come  up  with  a  guaranteed  optimum  in  the  allotted  time  -  in  feet 
all  it  could  guarantee  waa  that  the  solution  found  was  within  6.67%  of  optimum. 
After  more  than  double  our  normal  resource  maximum  was  expended,  and  the  MIP  had 
gone  through  306  iterations,  no  better  solution  was  found.  All  we  learned  was 
that  the  solution  already  found  was  guaranteed  to  be  within  6.48%  rather  thar 


6.67 %  of  optimum.  An  additional  MIP  run,  witn  a  different  APEX  option  (different 
search  procedure)  did  find  a  solution  which  reduced  the  MIP  objective  function 
by  0.7%.  In  addition  to  the  work  reported  here,  the  radio  had  been  previously 
run  a  number  of  times  as  the  MIP  FORKULATOR  pregram  evolved,  without  diffxculty. 

A  user  reports  no  problems  encountered  in  the  analysis  of  a  successor  to  the 
tactical  radio,  albeit  one  with  fewer  pieces  of  test  equipment. 

The  comparison  cf  heuristic  and  looping  solutions  are  found  in  Table  C.2. 
Because  of  the  problem  with  the  radio,  results  were  obtained  for  the  radio  with 
an  operational  availability  (OA)  target  of  95%,  a*  well  ae  for  the  target  of  90%; 
the  additional  target  gave  ua  another  test  case  to  investigate.  Only  in  the  case 
of  the  complex  radar  did  the  looping  solution  Improve  on  the  heuristic,  and 
then  only  in  terms  of  marginal  Improvement  In  availability  per  dollar. 

We  would  acknowledge  that  these  results  may  be  less  a  tribute  to  the  power 


of  the  heuristic  than  to  Insensitivity  of  the  maintenance  policies  chosen  to  the  \ 
used.  The  X  of  $0.01  for  the  simple  radar  Is  not  as  aberrant  as  it  may  seem,  since 
field  stockage  rules,  reflected  in  the  SESAME  atockage  routines,  dictate  that 
certain  levels  of  stockage  be  calculated  regardless  of  how  low  X  Is. 

One  other  surprising  aspect  of  the  results  vas  that  in  the  loop  for  the 
complex  radar,  the  X  values  kept  Increasing.  Intuitively,  we  had  expected  that  as 
we  increased  the  X  upon  which  the  MIP  inputs  were  based,  the  MIP  would  choose 
maintenance  policies  resulting  in  fewer  backorders,  so  that  the  X  required  in 
the  EVALUATOR  to  achieve  the  OA  target  would  decrease.  Under  the  correct  assump¬ 
tions,  this  intuition  can  be  buttressed  by  Lagranglan  theory. 

Here  is  what  caused  the  anomaly  for  the  complex  radar.  One  component  had  a 
maintenance  policy  change  during  the  looping  process.  At  the  lowest  X,  maintenance 
policy  (A)  projected  more  component  backorders,  but  was  cheaper  (including  backorder 
cost)  then  the  alternative  policy  (B) .  The  evaluator  raised  the  X  in  order  for 
policy  A  to  achieve  target  availability  by  reducing  the  component  backorders.  At 
this  next  higher  X,  the  MIP  chose  policy  B  over  A,  since,  even  though  policy  A  had 
fewer  component  backorders,  its  module  stockage  cost  was  high. 

Finally,  the  X  had  to  be  further  rained  in  the  evaluator  to  reduce  backorders 
of  the  now  chosen  policy  A  (which,  by  the  way,  was  a  mixed  policy;  hence  the  pro¬ 
ration  scheme  in  Section  3.4  was  used  in  the  MIP).  At  this  last  X  (31,059,002), 
a  final  pass  thru  MIP  and  evaluator  stabilised  on  policy  £  at  a  cost  of  $9,598,000. 
Note  that  this  .2%  Increase  in  coat  improved  the  unavailability  by  .6%  (from  .0148 
to  .0147). 


TABLE  C.l 


EXAMPLE  CHARACTERISTICS 
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EXAMPLE 

COMPONENTS 

NUMBER  OF 

MODULES 

FAILURE 

MODES 

TACTICAL  RADIO 

19 

61 

61 

SIMPLE  RADAR 

3 

16 

16 

COMPLEX  RADAR 

11 

37 

37 

$31,059,002  $31,059,002  $9598  98.53Z 


APPENDIX  D 


COMPUTER  COSTS  AND  RUN  TIMES 

Before  discussing  run  times  and  costs,  we  will  devote  a  section  to  a 
theoretical  principle  which  avoided  combinatorial  problems  In  programming 
the  assessment  of  maintenance  policies  over  many  Items. 


Principle  of  Independence  In  Use  of  SESAME  Subroutine 


One  principle,  that  contained  In  the  SESAME  optimization  process,  in¬ 
herently  reduces  the  run  times  In  the  MIP.  On  the  surface  one  might  think 
that  in  any  one  cost  evaluation  step,  specific  maintenance  policies  on  each 
LRU  would  have  to  be  considered  simultaneously  In  order  to  assess  cost  and 
performance  for  a  specified  Lagranglan  value;  therefore  the  number  of  com¬ 
binations  of  policies  over  the  LRUs  would  be  enormous. 

However,  In  the  SESAME  optimisation  process,  given  the  Lagranglan  value, 
the  Individual  item  supply  availabilities  are  "set"  based  on  the  Items'  unit 
prices  relative  to  this  value.  Each  item  (LRU)  la  optimized  Independently 
(multiechelon  allocation  of  stockage  quantities)  to  achieve  its  target. 

An  Item's  optimal  stockage  of  course  will  be  dependent  on  its  maintenance 
policy,  but  Independent  of  the  other  Items'  policies.  This  principle  allows, 
for  a  set  Lagranglan  value,  a  SESAME  subroutine  of  the  ORLA  program  to 
compute  the  optimal  kage  costs  for  an  Item  under  some  pure  maintenance 
policy  (e.g.  100Z  replacement  at  some  echelon  and  100X  repair  at  some  echelon) 
without  having  to  simultaneously  consider  hybrid  combinations  of  policies 
across  all  items.  The  ORLA  "optimizer"  uses  the  results  of  the  SESAME  module 
(stockage  costs,  backorders  by  Item  by  pure  policies),  in  conjunction  with 
other  logistical  costs,  to  find  the  best  maintenance  policy  vector  across 
all  items  for  a  given  Lagranp* . 


>erience  With  Com^uwer  Cc.'  .a 


A  commercial  time  sharing  service  was  used  for  the  work  reported  here. 

Costs  reported  are  based  on  privity  02  which  guarantees  12  hours  response  time. 
Actual  response  times  are  most  .an  a  few  minutes  and  never  exceed  a  few  hours. 
The  costs  reflect  a  substantial  volume  discount.  They  do  not  include  an  input/ 
output  charge  as  we  did  not  keep  an  accurate  record  of  these.  For  the  MIP  runs, 
costs  Include  a  large  surcharge  for  use  of  the  APEX  program;  this  surcharge 
is  computed  as  a  multiple  of  usage  costs. 
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The  cost  of  running  the  MIP  FORMULATOR  for  the  radio  was  about  $6.50.  Costs 
for  the  EVALUATOR  were  smaller.  This  was  in  line  with  expectations  since  the 
SESAME  routines  which  are  the  heart  of  these  programs  run  fast. 

The  range  of  costs  encountered  In  running  MIP  for  the  radio  is  shown  below. 
There  were  97  constraints  and  938  variables  of  which  30  were  integer.  Of  the 
30  Integer  17  were  binary,  and  the  largest  upper  bound  for  any  Integer  variable 
was  62.  The  different  runs  varied  only  in  the  coefficients  of  the  objective 
function.  The  sensitivity  of  run  costs  to  changes  in  coefficients  has  been 
noted  by  many  other  KIP  users. 


RUN 

COST 

Heuristic  -  OA  ■ 

90* 

$38 

Loop  -  OA  *  90X* 

$42 

Heuristic  -  OA  * 

95* 

$13 

Loop  -  OA  ■  95* 

$13 

Reducing  Running  Times  for  the  MIP 


The  logistician  often  knows  a  priori  that  certain  solutions  would  be  too 
costly  or  would  otherwise  be  unacceptable.  He  can  input  this  knowledge  in  two 
ways:  he  can  specify  the  lowest  echelon  at  which  a  piece  of  test  equipment  can 
be  placed,  or  he  can  rule  out  certain  policies;  e.g.  no  end  item  repair  above 
user  level. 

Preliminary  analysis  may  be  used  to  identify  upper  bounds  on  echelon  test 
equipment  requirements;  i.e.  the  requirement  if  all  repair  using  that  equipment 
were  performed  at  that  echelon.  Typically,  the  bound  will  be  1,  transforming 
the  integer  into  a  binary  variable. 

This  concept  can  be  carried  further  with  some  degradation  in  the  optimality 
of  the  solution.  Por  each  t^,  associate  t*  where  t'^  is  continuous  and 
let  M  be  a  "large"  integer. 

Set 

^ek  ±  (tekJ  3 


(M) 


or 


*'*  Mj.i  00 


Run  stopped  because  limit  on  computer  costs  reached. 


In  the  first  case  we  permit  continuous  values  for  test  equipment  provided  they 
exceed  1.  In  the  second  cess  continuous  values  less  then  1  ere  permitted 
If  et  least  1  teat  equipment  per  facility  has  been  deployed  at  a  lower  echelon. 
The  rationale  is  that  if  we  must  round  up  to  1  at  a  higher  echelon,  we  have  not 
miscalculated  too  severely  since  the  algorithm  alresdy  had  costed  out  Integer 
quantities  at  a  lower  echelon  where  there  are  typically  winy  more  facilities. 
Furthermore,  an  optimum  policy  would  not  tend  to  purchase  test  equipment  at  a 
higher  echelon  when  all  requirements  could  be  met  by  equipment  already  purchased 
for  the  lower  echelon. 

Under  the  simplified  formulation  we  may  add  the  constraint 


t  .  <  1 

ek  — 


which  further  speeds  the  MIP.  This  is  not  an  equality  because  of  the  possibility 
of  component  or  module  discard  in  which  case  the  equipment  is  not  needed  at  all. 
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